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^S^ ' We present a new analysis of the AI = 1/2 rule in K — > tttt decays and the 

Bk parameter. We use the 1/N C expansion within the effective chiral lagrangian for 
pseudoscalar mesons and compute the hadronic matrix elements at leading and next- 

*^ to-leading order in the chiral and the 1/N C expansions. Numerically, our calculation 

reproduces the dominant AI = 1/2 K — ► tttt amplitude. Our result depends only 
moderately on the choice of the cutoff scale in the chiral loops. The AI = 3/2 
amplitude emerges sufficiently suppressed but shows a significant dependence on the 
cutoff. The Bk parameter turns out to be smaller than the value previously obtained 
in the 1/N C approach. It also shows a significant dependence on the choice of the 
cutoff scale. Our results indicate that corrections from higher order terms and/or 
higher resonances are large for the AI = 3/2 K — * tttt amplitude and the (|A,S| = 2) 
K° — K° transition amplitude. 
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1 Introduction 

Over the last few decades the kaon system has provided us with a rich field of phenomenol- 
ogy which has been important for developing our theoretical understanding of the interplay 
of weak and strong interactions. The nonleptonic kaon decays are especially interesting be- 
cause they provide a testing ground for QCD dynamics at long distances. Two outstanding 
problems in the field are the explanation of the AI = 1/2 rule in K — > tttc decays and the 
calculation of the Bk parameter which measures the non-perturbative contributions to the 
{\AS\ = 2) K° — K° transition amplitude. An accurate knowledge of Bk is necessary for 
theoretically investigating the indirect CP violation in the neutral kaon mass matrix, as 
well as, the Kl — K$ mass difference. The AI = 1/2 rule is particularly important because 
it gives rise to the small value of the ratio e'/e which measures the direct CP violation in 
the K — > 7T7T decay amplitudes. 

Since its first observation more than 40 years ago Jl| the AI = 1/2 enhancement 
has attracted a great deal of theoretical interest trying to find the dynamical mechanism 
behind the approximate isospin selection rule, in particular within the standard model. 
Experimentally, the ratio of the AI = 1/2 and AI = 3/2 amplitudes in K — > 7r7r decays 
corresponding to I = and I = 2 in the final state, respectively, was measured to be 

1 Reoo MK _> { „ n) ) = 

u Rea 2 Re(K -► (7rvr) /=2 ) v ; 

with Aj = aj exp(i8i) and Si the final state interaction phases. This result was particularly 
enigmatic before the advent of QCD when only the current- current operator Q 2 arising 
from the W exchange was included in the analysis and, consequently, Rea /Rea 2 was 
expected to be around one. With the establishment of QCD our understanding of the 
AI = 1/2 selection rule improved considerably. Using the operator product expansion, the 
K — > 7T7T amplitudes are obtained from the effective low-energy hamiltonian for \AS\ = 1 
transitions [2-4], 

K f s r - % e« E *(A0ft(A0 (A* < me) , (2) 

Ci{fi) = z t (n) + ryii^i) , r = -&/£„ , £ q = V* s V qd . (3) 

The arbitrary renormalization scale /i separates short- and long-distance contributions to 
the decay amplitudes. The Wilson coefficient functions Ci(fi) contain all the information on 
heavy- mass scales. For CP conserving processes only the z% are numerically relevant. The 
coefficient functions can be calculated for a scale jjl > 1 GeV using perturbative renormal- 
ization group techniques. They were computed in an extensive next-to-leading logarithm 
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4 5Z ^ M 9l QLl^d L , 
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Qe = 


-8 2j s L q R q R d L , 
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Qs = 
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-8 ^2 2 6q ^ L<lR V RdL ' 
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(7) 



analysis by two groups 0, ||. The local four-quark operators Qi(fj) can be written, after 
Fierz reordering, in terms of color singlet quark bilinears: 

Qi = 4 s L 7^ L ul'J^ul , 
<?3 = 4 ^ s L 7 M ci L gz^gL , 

Q5 = 4 ^ s L ^d L qR'Jf.qR , 

9 

Qi = 4 ^ -e g s L 7 M d L g B 7 M g B , 
9 

where the sum goes over the light flavors (q = u, d, s) and 

Qr,l = ^(l± 75 )g, e q = (2/3, -1/3, -1/3). (8) 

Qz, . . . , Qq arise from QCD penguin diagrams involving a virtual W and a c or t quark, 
with gluons connecting the virtual heavy quark to light quarks. They transform as (8l, 1r) 
under SU(3)l x SU(3) r and solely contribute to AJ = 1/2 transitions. Q 7 and Q 8 are 
electroweak penguin operators 0, [8] which are less important for the AI = 1/2 rule. Long- 
distance contributions to the amplitudes Aj are contained in the hadronic matrix elements 
of the four-quark operators, 

(Q^))! = (tttt, /| Qiifj) \K°) , (9) 

which are related to the 7r + 7r~ and 7r°7r° final states through the isospin decomposition 

(Qi)o = -^=(2(n + n-\Q l \K ) + (n°n \Q l \K )) , (10) 

V6 

(Qi} 2 = ^={(n+n-\Q t \K°}-(n\°\Q t \K }) = \J\^ + A Q t \K+) . (11) 

They are difficult to calculate but can be estimated using non-perturbative techniques 
generally for \x around a scale of 1 GeV. 

Major progress in the understanding of the AI = 1/2 rule was made when it was 
observed that the short- distance (quark) evolution, which is represented by the Wilson 
coefficient functions in the effective hamiltonian of Eq. ([2J), leads to both an enhancement 
of the 1 = and a suppression of the 1 = 2 final state. The octet enhancement || in 
the (Qi, Q 2 ) sector is dominated by the increase of z 2 when \i evolves from M w down to 



\x ~ 1 GeV, whereas the suppression of the AI = 3/2 transition results from a partial 
cancellation between the contributions from the Q\ and Q 2 operators owing to a destruc- 
tive Pauli interference in the K + — ► 7r + 7r° amplitude. Another important short- distance 
enhancement was found to arise in the sector of the QCD penguin operators, in particular 
for Ze, through the proper inclusion of the threshold effects (and the associated incomplete 
GIM cancellation above the charm quark mass) . Nevertheless, it was concluded that the 
perturbative QCD effects are far from sufficient to describe the AI = 1/2 rule and QCD 
dynamics at low energies must be addressed. The long-distance enhancement of the matrix 
elements of the QCD penguin operators over the matrix elements of Q± and Q 2 was first 
conjectured and estimated in Ref. || in the vacuum saturation approximation (VSA) [J10| . 
The VSA approach, however, fails completely in explaining the AI = 1/2 rule, and a more 
refined method for the calculation of the hadronic matrix elements is certainly needed. 

Due to the non-perturbative nature of the long-distance contribution, a large variety of 
techniques has been proposed to estimate it (for some recent publications see Refs. [11-16]). 
Among the analytical methods, the 1/N C expansion |17| (N c being the number of colors) 



associated with the effective chiral lagrangian is particularly interesting. In this approach, 
QCD dynamics at low energies is represented by the 'meson evolution' of the operators, 
from zero momentum to /1, in terms of the chiral loop corrections to the matrix elements 
||, 18]. The authors of Ref. [18] calculated the loop corrections to the matrix elements 



of Qi and Q 2 and included the gluon penguin operator Q 6 at the tree level, consistent 
with the 1/N C expansion. They obtained an additional enhancement and suppression of 
the AI = 1/2 and A J = 3/2 amplitudes, respectively, systematically continuing the octet 
enhancement in the (Qi,Q2) sector to the long-distance domain. Numerically, a 2 was 
reproduced with an accuracy of 70 to approximately 100 %, whereas a [for A QCD = 300 MeV 
and m s (l GeV) = 125-175 MeV] was found to be around 65 - 80 % of the measured value, 
suggesting that the bulk of the physics behind the AI = 1/2 rule in kaon decays is now 
understood. One might note that the agreement with experiment is not improved by 



including the next-to- leading order values for the Z{ [19| . 

In this article we present a new calculation of the hadronic matrix elements in K — > tttt 
decays in the 1/N C expansion for pseudoscalar mesons. The paper contains several im- 



provements over the original approach of Ref. [yj| which are conceptually and numerically 
important. One improvement concerns the matching of short- and long-distance contribu- 
tions to the amplitudes, by adopting a modified identification of virtual momenta in the 
integrals of the chiral loops. To be explicit, we consider the two currents or densities in the 
chiral representation of the operators to be connected to each other through the exchange 
of an effective color singlet boson, and identify its momentum with the loop integration 



variable. The effect of this procedure is to modify the loop integrals, which introduces no- 
ticeable effects in the final results. More important it provides an unambiguous matching 
of the 1/N C expansion in terms of mesons to the QCD expansion in terms of quarks and 
gluons. The approach followed here leads to an explicit classification of the diagrams into 
factorizable and non- f actor izable. Factorizable loop diagrams refer to the strong sector of 
the theory and give corrections whose scale dependence is absorbed in the renormalization 
of the chiral effective lagrangian. The non-factorizable loop diagrams have to be matched 
to the Wilson coefficients and should cancel scale dependences which arise from the short- 
distance expansion. In a recent publication together with W.A. Bardeen and E.A. Paschos 
20fl we used this method to calculate the hadronic matrix elements of Qq and Q$ which 



dominate the ratio e'/e. In this paper we focus on the CP conserving amplitudes which, 
to a large extent, are governed by the current-current operators Q± and Q2- 



In Ref. [18[| a mass scale replacing the complete dependence of the exact expressions 
on the meson masses was introduced in the chiral logarithms. Another improvement of 
this paper is that we investigate the exact expressions for the matrix elements using the 
matching prescription discussed above, i.e., we evaluate the complete finite terms from the 
non-factorizable diagrams. Moreover, we calculate the whole of the matrix elements, that 
is to say, we also take into account the subleading penguin operators. For consistency with 
Ref. |20] we also include the small effects of the singlet r] . In the numerical analysis we take 



special care to separate the different contributions. In particular, we discuss the effect of the 
final state interaction phases which were not taken into account in Ref. fl8| . Uncertainties 



arising from the short-distance part of the calculation are estimated by comparing the 
amplitudes obtained from the LO and the NLO Wilson coefficients, respectively. Finally, 
we also investigate the size of higher order corrections to the hadronic matrix elements to 
critically examine the stability of our results within the pseudoscalar approximation. 

In the second part of this work we investigate the matrix element of the (|AS| = 2) 
K° — K° amplitude in the 1/N C expansion following the same lines of thought. The 
introduction to this calculation we postpone to the beginning of Section 5. Our results for 
the K — > inr matrix elements were already discussed in part in Refs. PT] , E2J. For a more 



detailed presentation of the general method we refer the reader to Refs. |20| , |23 . 

The paper is organized as follows. In Section 2 we review the general framework of the 
effective low-energy calculation and discuss the matching of short- and long-distance con- 
tributions to the decay amplitudes. Then, in Section 3 we investigate the K — > tttc matrix 
elements. We show explicitly that the scale dependence of the factorizable loop diagrams 
is absorbed in the renormalization of the bare couplings, the meson wave functions and 
masses. We next calculate the non-factorizable loop corrections in the cutoff regulariza- 



tion scheme. In Section 4 we match them to the Wilson coefficients to obtain the isospin 
amplitudes. In Section 5 we extend the analysis to the (| AS\ = 2) K° — K° transition. We 
compute the matrix element and match it to the short-distance coefficient function to de- 
termine the Bk parameter. In both sections we present our numerical results and compare 
them with those of the existing analyses. The conclusions can be found in Section 6. 

2 General Framework 



Following the lines of Ref. [2(J we calculate the hadronic matrix elements of the local four- 
quark operators (with \AS\ = 1, 2) in the 1/N C expansion. To this end we start from the 
chiral effective lagrangian for pseudoscalar mesons which involves an expansion in momenta 
where terms up to 0(p A ) are included p4[ . Keeping only terms of 0(p 4 ) which contribute 
to the K — > tttt or the K° — K° matrix elements and are leading in N c it reads:[] 



C 



eff 



P 

4\ l " ' AN C 

+rL 5 (D ll U j D^U(M j U + U ] M)) + r 2 L s (M ] UM ] U + MU ] MU ] ) , (12) 



((D^D^U) + ^-(InU^ -\nU) 2 + r(MU j + UM j ) 



with D^U = d^U — irJJ + iUl^, (A) denoting the trace of A and M. = diag(m u , m d , m s ). 
l^ and r M are left- and right-handed gauge fields, respectively, / and r are free parameters 
related to the pion decay constant F n and to the quark condensate, with r = —2(qq)/f 2 . 
The complex matrix U is a non-linear representation of the pseudoscalar meson nonet: 







U = exp - 


n, n = vr a A a , 
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in terms of the physical 


states 






f n° + -±ar] + - 


f\brf V2n+ 
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V2n- 


-7T + ±a V + 






{ V2K- 


V2K° 


with 












a = cos 9 — 


V2sm9, V2b = 



(A a A 6 ) = 28, 



ab 



(13) 



hi] 1 



\ 



-^+V3 a? ?' / 



V / 2 J ff+ 
V2K° 



(14) 



sin 9 + v2 cos 9 , 



(15) 



The various conventions and definitions we use are in agreement with Ref. |2(|. In partic- 
ular, we introduce the singlet r] in the same way and with the same value for the Ua{X) 

1 One might note that the mass term ex Lg contributes only to the matrix elements of Qe and Qs which 
were computed in Ref. pQj. Here we include it for completeness. 



symmetry breaking parameter, a = m 2 + mr, — 2m K ~ 0.72 GeV , corresponding to the 
T] — rf mixing angle 9 = —19° |[25| . The bosonic representation of the quark currents is 
defined in terms of (functional) derivatives of the chiral action: 

re £2 

+irL s (& l lflM - M j d^U + d^UM^U - tfMU*d' t U) ji , (16) 



and the right-handed currents are obtained by parity transformation. Eq. fll6"|) allows us 
to express the current-current operators in terms of the pseudoscalar meson fields. 

The 1/N C corrections to the matrix elements (Qi)i are calculated by chiral loop dia- 
grams in line with Ref. p0| . The factorizable contributions, on the one hand, refer to the 
strong sector of the theory and give corrections whose scale dependence is absorbed in 
the renormalization of the chiral effective lagrangian. This property is obvious in the case 
of the (conserved) currents and was demonstrated explicitly in the case of the bosonized 
densities [^0[ |^j. Consequently, the factorizable loop corrections can be computed within 
dimensional regularization. The non-factorizable corrections, on the other hand, are UV 
divergent and must be matched to the short- distance part. They are regularized by a fi- 



nite cutoff which is identified with the short-distance renormalization scale 0, [19], [26], £7| . 
The definition of the momenta in the loop diagrams which are not momentum translation 
invariant was discussed in detail in Ref. ||20|| . A consistent matching is obtained by con- 



sidering the two currents or densities to be connected to each other through the exchange 
of a color singlet boson and by assigning the same momentum to it at long and short 
distances [28-31]. The identification of this momentum with the loop integration variable 



leads to modified integrals in the chiral loop diagrams compared to those of Refs. [18], ^6| . 
The numerical implications for the isospin amplitudes in K — ► mr decays and the Bk 
parameter will be addressed in Sections 4 and 5. 

In this paper we investigate the hadronic matrix elements at leading and next-to-leading 
order in the chiral and the 1/N C expansions. In particular, we calculate the 0(p 2 /N c ) 
corrections to the current- current operators, that is to say, the one-loop corrections over 
the 0(p 2 ) lagrangian. The matrix elements of the density-density operators Q§ and Qg 



are taken from Ref. [ 2"0"[ . In the numerical analysis of the A/ = 1/2 rule and the B^ 
parameter we use the leading logarithmic (LO), as well as, the next-to-leading logarithmic 
(NLO) values || || |32], |33|] for the (\AS\ = 1, 2) short- distance coefficient functions.^ In 



2 We treat the coefficient functions as leading order in 1/N C since the large logarithms arising from 
the long renormalization group evolution from (mt,Mw) to fi ~ 0(1 GeV) compensate for the 1/N C 
suppression. 



general, the lack of any reference to the renormalization scheme dependence in the effective 
low-energy calculation prevents a complete matching at the next-to-leading order [EI . 



Nevertheless, a comparison of the amplitudes obtained from the LO and NLO coefficients 
is meaningful in order to test the validity of perturbation theory. 

In the following sections we calculate the long-distance 1/N C corrections to the K — > txtx 
amplitudes and the Bk parameter. First, we investigate the factorizable corrections and 
show their absorption in the low-energy constants. Secondly, we determine the non- 
factorizable loops within the modified momentum prescription. Finally, we perform a 
numerical analysis and compare our results with those of the existing studies. 

3 K — > 7T7T Decays 

In this section we present the hadronic matrix elements of the current-current operators for 
the physical decay modes K° — > tt + tt~ and K° — > 7r°7r° up to 0(p A ) and 0(p 2 /N c ) in the 
parameter expansion. From these results we derive the isospin amplitudes K — > (7T7r)/ = o,2, 
heading for an explanation of the AI = 1/2 selection rule in kaon decays. 

3.1 Factorizable 1/N C Corrections 

The (bare) tree level of the K — > titi matrix elements, up to 0(p 4 ) in the chiral expansion, 
as well as, the factorizable 1/N C corrections to the 0(p 2 ) can be calculated from the tree 
and loop topologies depicted in Fig. [I]. From the sum of these diagrams we obtainP] 



(n+n-\Q 2 \K% } = V2 f (m 2 K - ml) 



4L 5 



1 + -^ (m 2 K + 4ml) 



16w 2 f 2 
where 



3A^K + 2m;;)logA c 2 ) + 



(17) 





(n + n-\Q 2 \K°) F = (n + n-\Q 4 \K°) F 


= -(vrVlQ^y 




= (n°n°\Q,\KY 


9 

= g (tt 71 \Q7\X ) 


and 








(n+n-lQilKY = for 


i G {1,3,5,7} 




(tt° tt° \Qi\K°) F = for 


% E {2,3,5}. 



(19) 

(20) 

3 In distinction to Ref. |2(J the factor i referring to the weak vertex is included in the definition of the 
matrix element. 
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(ij)(kl) 





(ij)(kl) 








mm 



S^ (ij)(kl) , (ij)(kl) 



+ (ij)^-^(kl) 



Figure 1: Factorizable diagrams for the matrix elements of the current-current operators 
in the isospin limit. Crossed circles represent the bosonized currents, black circles the 
strong vertices. The lines denote the pseudoscalar mesons. The external legs represent all 
possible permutations of the kaon and the pions. 



The ellipses in Eq. (|17|) denote finite terms we omit here for the analysis of the ultraviolet 
behaviour (in particular, they provide the reference scale for the logarithms). We spec- 
ify our results in the cutoff regularization scheme to demonstrate the absorption of the 
quadratic, as well as, the logarithmic divergences as required by current conservation. We 
note that all factorizable terms quadratic and logarithmic in the cutoff are independent of 
the momentum prescription in the loop. A c is the cutoff for the factorizable diagrams. We 
introduce two different scales since the factorizable and the non-factorizable corrections re- 
fer to disconnected sectors of the theory (strong and weak sectors). Having demonstrated 
the absence of UV divergent terms in the sum of the factorizable diagrams, in the numeri- 
cal analysis of the full expressions we will use dimensional regularization, as in pure chiral 

perturbation theory, which is momentum translation invariant. 

1 ii 
If we renormalize the wave functions of the kaon and the pions (7r r = Zj tto), as well 

as, the bare decay constant / by using Eqs. (14)-(17) and (25) of Ref. p0| , we arrive at the 

renormalized (factorizable) matrix elements of the (\AS\ = 1) current-current operators:^ 



(n+n-\Q 2 \K°) F 



(r) 



V2F W (m 



K 



ml) 



1 + 



4L| 
F 2 



m„ 



(21) 



where the constant L\ is defined through the relation [20 



F. ~ Fl 



nr 



K 



m„ 



(22) 



4 The full expressions for the wave function and the decay constants are given in terms of integrals in 



Appendix A of Ref. |20 






(ij)(kl) ^ (ij)(kl) y (ij)(kl) 



+ (ij)-«-Mkl) 



Figure 2: Non-factorizable diagrams for the matrix elements of the current- current oper- 
ators in the isospin limit. 




and the remaining matrix elements can be obtained from Eqs. fll8D-([20|). 

We notice that for the four-quark operators Qi of the current- current type the divergent 
terms are absorbed by the renormalization procedure. In addition, the factorizable 1/N C 
corrections vanish completely, that is to say, the divergent as well as the finite terms. 
This property has been observed numerically, within dimensional regularization, because 
the complexity of all factorizable contributions prevents us from doing a fully analytic 
calculation. Since the factorizable scale A c disappears through renormalization, the only 
matching between long- and short- distance contributions is obtained by identifying the 
cutoff scale A c of the non-factorizable diagrams with the QCD renormalization scale. 

Finally, we note that in the next-to-leading order term of Eqs. ( plf ) and (|22"D we used 



1/F-n rather than 1/ / as it was done in Ref. []T8[] . Formally, the difference represents higher 
order effects. Nevertheless, the appearance of 1/ / gives rise to a residual dependence on 
the factorizable scale A c , which has no counterpart at the short- distance level and will be 
absorbed by factorizable loop corrections to the matrix elements at the next order in the 
parameter expansion. Consequently, it is a more adequate choice to use the physical decay 
constant in the expressions under consideration. Instead of F n the kaon decay constant 
Fk could be used as well. Both choices will be considered in the numerical analysis, which 
gives a rough estimate of higher order corrections. 

3.2 Non-factorizable 1/N C Corrections 

The non-factorizable 1/N C corrections to the hadronic matrix elements constitute the part 
to be matched to the short-distance Wilson coefficient functions; i.e., the corresponding 
scale A c has to be identified with the renormalization scale \x of QCD. We perform this 



10 



identification, as we argued in Section 2, by associating the cutoff to the effective color 
singlet boson. Then, at the 0(p 2 ) in the chiral expansion, from the diagrams of Fig. ^| we 
obtain in the 577(2) limit: 



(n + n-\Q 2 \K°) NF 
(n+n-\Q 3 \K°) r ~ 

(n + n-\Q A \K°) NF 

(n + n-\Q 7 \K°) NF 



y/2 (m 2 K - m\ 



16tt 2 F t 


y/2 {m 2 K - m 2 n ) 


16vr 2 F T 


\/2 {m\ - ml) Q 


16vr 2 F T 


y/2 (m\ - ml) 


1Q7T 2 F V 


V2(m 2 K + 2ml) 



3 A? 



-m 2 K + 3ml) logA^ + 



^A 2 c+ (m 2 K -^ml) leg A; 



2mi logA^ + -- 



9 
2 A ' 



5 
2 



A?+f7<-^m!] logA 2 ^ 



167r 2 77 



6 m, 



A 2 C -- (3m 2 K + 7ml+ , , ., 
I 8 V m|. + 2 m 2 



logA 2 + -- 



0\JVF 



K7T U |Q 2 |K U > 

(7r 7r |g 4 |K°) 7VF 
where 



V^K,- 


m 2 ) 


16vr 2 F^ 




V2 (m 2 K - 


m 2 ) 



16tt 2 F^ 



2 A ' + 4 (m^-6rr4) logA c 2 + 



^A^ + ^Sm^-lOm^logA 2 



(7r + 7r-|g 3 |^ U ) 



Q\NF 



n»n»\Q 3 \K») 



0\NF 



1 



-(7T 7T |g 7 |K U ) 



Q\NF 



[n°n \Q,\K°) NF 
-(n+n-\Q 5 \K°) NF 



and 



(TrVlQil/T')^ = 0. 



(23) 
(24) 
(25) 
(26) 



(27) 
(28) 
(29) 



(30) 



(31) 



One might note that in Eqs. (|23])-(|2P|) [as in Eq. flTT|)] we replaced m 2 m 2 ,, and the mixing 
angle 9 by ml and m 2 K using the octet-singlet mass matrix of Ref. | f25f . 

At this stage of the calculation we find quadratic, as well as, logarithmic divergences 
of the non-factorizable corrections. We note that already the leading (~ A 2 ) terms depend 
on the momentum prescription. The quadratic terms were calculated in Ref. [3D|| in the 
background field method. In this paper we investigate the full expressions for the matrix 



11 



elements needed for the numerical analysis of the amplitudes. The results contain finite 
terms, originating from the solutions of the integrals listed in Eq. fl4"5|) and in Appendix B 
of Ref. |[20|| , which we neglect here for brevity and denote by the ellipses. We also note 
that in the case of Q? the solution of the integrals brings along a quartic dependence on 
the cutoff which has to be cancelled by adding a specific contact interaction proportional 
to S^'(0) to the Feynman rules of the truncated meson theory J30, |35fl . 



Even though the scale dependence of the perturbative coefficient functions is only loga- 
rithmic, the full long-distance contribution including the quadratic terms has to be matched 
to the short-distance part. The quadratic dependence on the cutoff is physical and is neces- 
sary for several reasons. First, in the chiral limit (m q = 0) all corrections vanish except for 
the A? c terms, which bring in the only scale to be matched to the short distance. Secondly, 
they stabilize the 1/N C expansion and generally improve the matching of the meson and 
the quark pictures |18| . Finally, they provide us with a rough estimate of the contributions 
from higher resonances. 

We note that in Eqs. (|23|) (p9|) we used the physical decay constant F n rather than / 
in the same way as for the factorizable diagrams. Again the difference represents higher 
order effects. However, the (factorizable) scale dependence of / has no counterpart in the 
short distance and will be absorbed at the next order in the chiral expansion. As for the 
factorizable contributions the choice of Fr instead of F n would be also appropriate. 

4 Numerical Analysis 

In this paragraph we list the numerical values for the hadronic matrix elements. We next 
match them to the Wilson coefficients and study the K — > (mv)i = o t 2 isospin amplitudes. 



In Section |4.1| we discuss in detail the 1/N C corrections to the matrix elements. In this 
context we also calculate the bag parameters, which quantify the deviations from the 
results obtained in the vacuum saturation approximation and, therefore, are convenient 
for a comparison with other works. The main results of the present analysis can be found 



in Section 4.2. Therein we give the amplitudes a® and 02 as functions of the matching scale 



and compare them with the data. 

4.1 Hadronic Matrix Elements 

Throughout the numerical analysis we use the following values for the parameters 

m n = (m n o + m 7r +)/2 = 137.3 MeV, F n = 92.4 MeV, 

12 



m K = (m K o + m K +)/2 
m v = 547.5 Me V, 
Try = 957.8 MeV , 



2 = 495.7 MeV, F x = 


113 MeV, 


e 


-19°, 


G F = 


1.1664- 1(T 5 GeV" 2 


\v us \ = 


0.22. 


we compute L r 5 = 2.07 x 10 -3 . 





\V ud \ = 0.974, 

Substituting them in Eq. (|2!_ 

We parameterize our results in terms of the non-perturbative bag parameters B\ and 

(3/2) 

B\ , which quantify the deviations from the values obtained in the vacuum saturation 
approximation fi~0|| : 



B, ,1/2) = HI?, *e{l,..,8>, (32) 

Bf m = WM, i 6 {1,2,7,8}, (33) 



'2 



with (Qi)i containing both factorizable and non-factorizable contributions. The VSA ex- 
pressions for the matrix elements are taken from Eqs. (XIX.11)-(XIX.28) of Ref. [fTT)|.p| 
The numerical values for the matrix elements of the current-current operators are given in 
Tables |l| and [| (Qs)^ SA and (Q?)^ A are functions of R = 2m 2 K /(m s + m d ) ~ 2m 2 K /m s 
and, consequently, depend on the renormalization scale. For comparison, in the tables we 
also show the results obtained in the large- N c limit, see Eqs. (|i~8|)-(f21"|). One might note 
that the different values generally do not coincide, even if the small 0{p A ) term propor- 
tional to m 2 in Eq. ([H]) [which contributes only at the level of 2 % of the 0(p 2 ) tree level 
term] is neglected, since in the vacuum saturation approximation Fierz terms are taken 
into account which are subleading in N c . In particular, the matrix element (Qi)^ SA differs 
by a factor of (1 — 2/N c ) from the result obtained at the 0(p 2 ) in the large-iV c limit. We 
notice that the inclusion in part of the 1/N C corrections in the VSA method leads to a 
suppression and enhancement of the 1 = and 1 = 2 amplitudes, respectively, in complete 
disagreement with the data. 

In Tables |3| and |4] we list our results for the hadronic matrix elements at next-to- 
leading order in the chiral and the 1/N C expansions. The matrix elements of the current- 
current operators are calculated from Eqs. (|i~8D-(|31|) including the finite terms denoted by 



the ellipses. The results for the operators Qq and Qs are taken from Ref. |p0| . These 
results contain the leading plus next-to-leading order terms in the chiral expansion of the 
density-density operators, as well as, the leading l/N c corrections, that is to say, the O(p ), 

5 Note that our definition of the pion decay constant (F^ = 92.4 MeV) differs by a factor of l/v2 from 
the one used in Ref. |19[ . 
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(Qi)o 


(Qa)o 


(Qsh 


(Qa)o 


(Q 5 )o 


(Qt}o 


VSA 
tree 


-4.03 
-12.3 


20.2 
24.6 


12.1 



36.3 
37.0 


-11.7 -R 2 



18.2 + 32.5 • R 2 
18.5 



Table 1: 1 = matrix elements of the current- current operators: VSA vs. tree level 
(large- N c limit), in units of 10 6 • MeV 3 (R in units of GeV). 





(Qih 


(Q 2 ) 2 


(Qrh 


VSA 
tree 


22.8 
17.4 


22.8 
17.4 


-25.7+ 18.9 -R 2 
-26.1 



Table 2: Same as in Table FLl now the 1 = 2 matrix elements. 



0(p 2 ), and O(p /N c ). Note that the matrix elements generally contain a non- vanishing 
imaginary part (scale independent at the one-loop level) which is due to on-shell (ir — it) 
rescattering effects. 

The isospin amplitudes are largely dominated by the operators Q± and Qi. Therefore 
it is instructive to analyze in detail the 1/N C corrections to these two operators. To this 
end we next give the analytic expressions for the isospin matrix elements of Q\ and Q 2 : 



(Qi)i 



l 



F n (m 2 K - ml) 



, 4L 5 2 

1 + -rSml + 



1 



n w w 2 *? 



x ( 6A^ 



l -m 2 K + Qml) log (l + ^ 



+ Qio \m\ 



(Q 2 )c 



— =F n [m K - mj 



-, 4L 5 2 

1 + -^W + 



F 2 * {An?F 2 



x I -A 2 + f "m| - ^ml) log (l + |§ 



15 
T 






(Qi> 2 = (Q 2 ) : 



i 4L 5 2 1 

1 + -Sml + 



x 



(-3A? +(!< + 3mj) log (! + §)) 



+ 021 m 



(34) 



a 20 [m] (35) 



(36) 



Eqs. (|34|)-(p6D allow us to compare our results with the analytic expressions of Ref. |L8 
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A c 


0.5 GeV 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 




<Qi)o 


-27.4 


-33.2 


-40.2 


-48.2 


-57.3 


-67.4 


-5.55i 


(Q 2 )o 


50.0 


58.8 


68.8 


79.9 


92.4 


106 


11. U 


(Q 3 )o 


0.04 


0.05 


0.03 


-0.02 


-0.12 


-0.26 





(g 4 ) 


77.5 


92.1 


109 


128 


150 


173 


16. 6i 


(g 5 } 


-0.04 


-0.05 


-0.03 


0.02 


0.12 


0.26 





(Q 6 )o 


-44.1 


-38.6 


-33.7 


-29.4 


-25.5 


-21.9 





(Q 7 )o 


34.4 


40.1 


46.6 


54.1 


62.6 


72.2 


8.32i 


(Q 8 )o 


118 


119 


119 


119 


118 


117 


36. 1% 



Table 3: Hadronic matrix elements of Qi,...,5,7 (in units of 10 6 • MeV 3 ) and (^6,8 (i n units 
of R 2 ■ MeV) in the isospin limit for the 1 = amplitudes, shown for various values of the 



cutoff A c . 



A c 


0.5 GeV 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 




(Qih 


6.54 


2.51 


-2.26 


-7.77 


-14.0 


-21.1 


-3.45? 


(Ga>2 


6.54 


2.51 


-2.26 


-7.77 


-14.0 


-21.1 


-3.45z 


(Qih 


-14.5 


-10.7 


-6.27 


-1.15 


4.67 


11.2 


5.18z 


(Qsh 


39.9 


35.3 


31.2 


27.2 


23.2 


18.8 


— ll.Si 



Table 4: Same as in Table EjL now for the 1 = 2 amplitudes. 



First, we note that the modified matching which was discussed in Section 2 increases the 
terms quadratic in the cutoff by a factor of 3/2 relative to the results presented therein. 



This was already observed in Ref. [30]. The modification of the quadratic terms provides 
an additional octet enhancement in the long-distance domain. The logarithmic terms, on 
the other hand, are modified only on account of the presence of the tjq. To be explicit, in 
the octet limit [i.e., in the absence of the i] , with a = b = 1 and m 2 = (Am 2 K — m 2 )/3] 
the coefficient of the logarithm in Eq. (|3"4f) is reduced to (m 2 K /2 + 10m^/3) whereas the 
other terms remain unchanged. The separation of the logarithmic and the finite terms in 
Eqs. (J34|)-( l36l) is arbitrary and is done, for comparison with Ref. |0|, by introducing a mass 
scale replacing the dependence of the exact expressions on the meson masses in the chiral 
logarithms. The logarithmic and the finite terms (an) defined in this way each depend 
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(Qi)o 


(Ql>2 


(Q 2 )o 


(Q 2 } 2 


tree 


-12.3 


17.4 


24.6 


17.4 


A? 


-34.5 


-24.4 


43.1 


-24.4 


logA c [m] 


4.43 


3.13 


10.0 


3.13 


finite 


-5.83 - 5.55i 


-3.90 - 3.45i 


2.20 + 11. li 


-3.90 - 3.45i 


total 


-48.2 - 5.55i 


-7.77-3.45i 


79.9 + 11. U 


-7.77-3.45z 



Table 5: Different contributions to the hadronic matrix elements of Q\ and Q2 (in units of 
10 6 • MeV 3 ) for A c = 800 MeV and m = 300 MeV. 

on the choice of the mass scale m, whereas the sum of all contributions is independent of 
this parameter. We calculated the complete finite terms arising from the non-factorizable 
loop diagrams using the matching prescription advocated in Refs. |20, 0.0 These terms 



were not included in Ref . |TJ| . Consequently, the numerical values of the matrix elements 
reported therein exhibit a dependence on the specific choice of the mass scale in the loga- 
rithms which is absent in the present calculation. 

In Table [5] we split up the numerical values for the 1 = and 1 = 1 matrix elements of 
Qi and Qi with respect to the quadratic, the logarithmic, and the finite terms, respectively, 
at a cutoff scale of A c = 800 MeV. From the table we see that the finite terms are of the 
same order of magnitude as the logarithmic ones and, therefore, must be considered at the 
same level in the numerical analysis. These terms are generally suppressed by a factor of 



S 



m 



K,TT 



/(A.'kFk) 2 < 20 % with respect to the leading 0{p 2 ) tree level. In addition, as can 



be seen from Eqs. (|~"~") -(""""") and Table ~5], no coefficient larger than one or two which could 
significantly enhance them has been found. This is different from the quadratic terms which 
are not suppressed as their relative size is determined by A = A 2 /(47rF n ) 2 and, moreover, 
they appear with larger prefactors [even as large as six in Eq. fl3~4Dl.fl Consequently, in 



6 For details on the computation of the loop integrals see Appendix B of Ref. ~~~|. 

7 It is interesting to note that the non-suppression of the quadratic terms presumably could be important 
for Qe but less important for Q%. On the one hand, the first non- vanishing tree level contribution to the 
operators Qe and Qs is of the 0(p 2 ) and O(p ), respectively. On the other hand, the first non- vanishing 
quadratic corrections to both operators are of the 0(p 2 /N c ) (terms of the O(p /N c ) were found to be 
only logarithmic fl2~|). Consequently, in the case of Qs the quadratic terms are (chirally) suppressed by a 
factor of p 2 ■ A with respect to the (leading) tree level contribution whereas in the case of Qq they bring 
in only a factor of A. Quadratic terms, even though subleading in N c , could therefore significantly affect 
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the case of the 1 = matrix elements of Qi and Q2 both the logarithmic and the finite 
corrections are moderate, and the chiral limit gives a satisfactory representation of the 
full amplitude provided that the matching scale is taken sufficiently large (A c > 500- 
600 Me V). In the case of the 1 = 2 matrix elements we also observe that the quadratic 
terms are enhanced with respect to the tree level, whereas the logarithmic and the finite 
terms are largely suppressed. However, in this case the quadratic corrections counteract the 
tree level, and the sum of both contributions is no longer large compared to the logarithmic 
and the finite terms. Therefore the neglect of either of the terms is no longer justified. 
In particular, we observe that for the A/ = 3/2 channel the chiral limit gives a better 
approximation to the exact result than a calculation which includes only the logarithms 
without taking into account the finite terms. This remark also holds for the matrix element 
(Qi)q. Finally, we note that variation of the mass scale in the logarithms [m^ < fh < ttlk} 
in Ref. fl8j has a noticeable effect on the numerical value of the 1 = 2 amplitude. 



When comparing the results of the present analysis with those of Ref. [[H| one has to 
take into account another difference in the treatment of the next-to-leading order terms: 
in Eqs. fl3"4|)-(|3~6"D we used 1/F n rather than the bare parameter 1/f as it was done in 



Ref. [0. Formally, the difference concerns higher order effects, as we already discussed 
above. However, since the factorizable scale which appears in the bare coupling / will 
be absorbed by factorizable loop corrections to the matrix elements at the next order in 
the parameter expansion, it has not to be matched to any short-distance contribution. 
Consequently, it is a more adequate choice to use the physical decay constant in the 
expressions under consideration. The effect of this different treatment of the next-to- 
leading order terms will be further discussed in Section PO . 



In Tables |6] and |7| we list the values we compute for the bag parameters B\ and 
B\ . We find a large enhancement of B x and B 2 over the VSA result, which 
constitutes the dominant contribution, at long distances, to the A/ = 1/2 transition in 
K — > 7T7T decays. Moreover, we obtain the correct scale dependence counteracting the 
scale behaviour of the Wilson coefficients z\ and Z2, which leads to an acceptable matching 
(see Section |4.2| ). In view of the large corrections one might question the convergence 
of the 1/N C expansion. However, there is no strong reason for such doubts because the 
non-factorizable contribution we consider in this paper represents the first term in a new 
type of a series absent in the large- N c limit. It is reasonable to assume that this leading 
non-factorizable term carries a large fraction of the whole contribution |Tj| (see also the 



the matrix element of Qe especially if large prefactors are observed as for Q\ and Q2 in Eqs. (|34|)-(36). 
This difference between the Q§ and Qs operators could play an important role for e'/e. This point will be 
investigated in Ref. |37] . 
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A c 


0.5 GeV 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


S (l/2) 


6.75 


8.24 


9.98 


12.0 


14.2 


16.6 


B am 


2.47 


2.91 


3.41 


3.96 


4.57 


5.23 


R (l/2) 


0.003 


0.004 


0.002 


-0.002 


-0.010 


-0.021 


d(1/2) 
• D 4 


2.12 


2.54 


3.00 


3.53 


4.13 


4.75 


B (l/2) 


0.0004 


0.0009 


0.0005 


-0.0003 


-0.0014 


-0.0020 


fl (l/2) 


1.26 


1.10 


0.96 


0.84 


0.72 


0.62 


B (l/2) 


0.15 


0.16 


0.18 


0.21 


0.23 


0.26 


o(l/2) 
-°8 


1.20 


1.21 


1.21 


1.21 


1.20 


1.19 



Table 6: Bag parameters for the 1 = amplitudes, shown for various values of the cutoff. 
B\ 7 § depend on R ~ 2m 2 K /m s and are calculated for a running m s (^ = A c ) at the leading 
logarithmic order (A QCD = 325 Me V) with m s (l GeV) = 175 MeV. 



A c 


0.5 GeV 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


^(3/2) 


0.29 


0.11 


-0.10 


-0.34 


-0.61 


-0.92 


^(3/2) 


0.29 


0.11 


-0.10 


-0.34 


-0.61 


-0.92 


£(3/2) 


-0.15 


-0.10 


-0.06 


-0.01 


0.04 


0.09 


R (3/2) 


0.72 


0.64 


0.56 


0.49 


0.42 


0.34 



Table 7: Same as in Table || now for the 1 = 2 amplitudes. 



? (l/2) 



|(l/2) 



discussion in Section f4~2| ). B^ and B§ turn out to be very close to zero. This property 
is due to the vanishing tree level, as well as, to the small 1/N C corrections proportional to 
171^/(4:71 F n ) 2 , see Eqs. (p5|) and (]3~0| ). We notice that the small contribution of the operator 
Qs to s'/e is even further reduced when replacing the VSA expression for (Qs)o, which is 



commonly used in the analysis of e'/e [Q, by the result presented in this paper. B 7 



(3/2*) 

and B\ are also found to be significantly reduced with respect to vacuum saturation 
approximation. In particular, B\ turns out to be negative for small values of the cutoff.Q 

8 Very recently [Q the first non-trivial \/N c corrections to the matrix elements of Q7 were evaluated 
using the methods of Ref. [p9| . The numerical results were also sensitive to the choice of the renormalization 



,(1/2) 



(3/2) 



scale. In particular, negative values for B\ and B) were found below /1 < 1.3 GeV, in qualitative 
agreement with the results of the present analysis but in disagreement with the large positive values 
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We also notice a decrease of the B\ and B\ parameters, which are relevant for A 2 . 
However, as we will see below, their scale dependence largely overcompensates for the 
variation of the short-distance coefficient functions. Nevertheless, as the values are found 
to be reduced, they generally account for the reduction of the 1 = 2 amplitude. Finally, 
Bq receives only small corrections whereas B% comes out to be substantially reduced 



relative to the VSA result [20 1. The numerical implications for e'/e will be investigated 
elsewhere [37] . One might note that the numerical values of Bg shown in Table [7| differ 



from the ones given in Table 2 of Ref . |2{J . This is due to the fact that in the present paper 



we include only the real part of the hadronic matrix elements in the definition of the Bi 



parameters (see Section [4~2l) . 



4.2 The AI = 1/2 Rule 

We next investigate the CP conserving amplitudes Rea and Rea 2 . To this end we start 
from the expression for the isospin amplitudes Ai which contain the (n — re) strong inter- 
action phase shift for the 1 = and the 1 = 2 final states, respectively, 

A 7=0 , 2 = ^ V ud V: s J2 otto (Q,(/i))/=o,2 ■ (37) 



Then 



Gr 



G 



Rea 7 = -| V ud V: s | £ z % (Q,)/ = -| V ud V* ua — £) * MQih ■ (38) 

* i * i 

Within an exact realization of non-perturbative QCD the two expressions in Eq. fl38"D are 
equivalent. However, in the approximate low-energy calculation of the present work the 
long-distance imaginary part which we computed at the one-loop level (see Tables |3| and f|) 
is not expected to be of the same accuracy as the real part obtained at this level. In par- 
ticular, as the one-loop (long-distance) imaginary part is scale independent, it cannot com- 
pensate for the scale dependence of the Wilson coefficients Zi leading to a scale dependent 
imaginary part of the total amplitude. This requires a calculation of the (long-distance) 
imaginary part at least at the two-loop level which will introduce a scale dependence. In 
addition, the two-loop contribution is expected to be of the same order of magnitude as the 
one-loop contribution which only appears at the level of the finite terms, as it will bring in 
a quadratically divergent term. This situation is analogous to the non- suppression of the 
one-loop contribution to the real part (~ A) with respect to the tree level. The two-loop 
contribution to the real part, on the other hand, is expected to be suppressed by at least a 

obtained in the chiral quark model at a matching scale of 0.8 GeV Eq |. 
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factor of S with respect to the tree level and the one-loop contribution. This is analogous 
to the one-loop logarithmic and finite terms which are suppressed by a factor of 5 with 
respect to the tree level. For the numerical analysis we will therefore consider only the real 
part of the matrix elements [see the second expression in Eq. (|38|) 1 using the experimental 
values of the final state interaction phases, 5 c xp = (37 ± 3)° and ^ xp = (-7 ± 1)° [|TJ. This 



procedure has also been followed in Ref. [42|. However, as the imaginary part is a loop ef- 



fect (suppressed by a factor of 5 with respect to the tree level contribution), its effect on the 
absolute value of the amplitude strictly speaking is of the two-loop order. Consequently, 
we will also compare our results with the ones obtained by taking the (long-distance) imag- 
inary part to zero, i.e., by taking J2i z i (Qi)i = Yl,i z i^ e (Qi)i- This holds for an estimate 
of the size of higher order effects which is generally disregarded in the literature. 

In Table |8| we show the numerical values of the amplitudes for various values of the 
matching scale and fixed values of A QCD = A^i- and the strange quark mass m s . The 
numerical analysis is done using the leading logarithmic, as well as, the next-to-leading 
logarithmic values of the Wilson coefficients listed in the appendix. The NLO values are 
scheme dependent and are calculated within naive dimensional regularization (NDR) and 
in the 't Hooft-Veltman scheme (HV), respectively^ The difference between the two NLO 
results at a given scale reveals the uncertainty due to the lack of any reference to the 
renormalization scheme dependence in the effective low-energy calculation. 

In Fig. |3] we show Rea calculated with leading order Wilson coefficients for various 
values of A QCD as a function of the matching scale. We take the (conservative) range of 
A QCD = 325±80 MeV which corresponds to a s (M z ) = 0.118±0.005 @. First, we note that 
our result for Oq shows an additional enhancement (around 30 - 50 % of the experimental 



value) compared to the result of Ref. JT8| which renders the amplitude in good agreement 
with the observed value for low values of the scale or even larger than the experimental 
value for large values of the scale. A significant enhancement arises from the Q± and 
Q2 operators due to the modified matching prescription in the non-factorizable sector we 
discussed above. Numerically, at a scale of A c = 800 MeV the modified momentum routing 
accounts for approximately 20 % of the final number(s) presented in Fig. 131. Another 



enhancement with respect to Ref. [Tj| originates from the correction of the real part by the 
experimental phase [see Eq. (J38"|)1. Neglecting completely the effect of the (it — ir) phase 
shift would reduce our result by a factor of cos S ~ 0.8. The remainder is due to the choice 
of the physical value F n instead of / in the next-to-leading order terms of the factorizable 



9 We are very thankful to M. Jamin for providing us with the numerical values of the Wilson coefficients 
used in this section. 
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Ac 


Reao 


Rea2 


LO 


NDR 


HV 


LO 


NDR 


HV 


0.5 GeV 


3.90 


0.74 


4.48 


0.063 


0.086 


0.063 


0.6 GeV 


3.50 


2.58 


3.57 


0.027 


0.032 


0.028 


0.7 GeV 


3.53 


2.89 


3.45 


-0.025 


-0.028 


-0.025 


0.8 GeV 


3.75 


3.13 


3.58 


-0.090 


-0.101 


-0.095 


0.9 GeV 


4.08 


3.42 


3.83 


-0.167 


-0.188 


-0.178 


1.0 GeV 


4.49 


3.76 


4.17 


-0.257 


-0.289 


-0.274 


exp. 


3.33 


0.15 



Table 8: Rea and Rea 2 (in units of 10" 4 MeV) for m s (l GeV) = 175 MeV, A. 
325 MeV, and various values of the matching scale \i = A c . 



A (4) 

QCD — ^ms 



and non-factorizable corrections. Our result depends only moderately on the matching 
scale although the stability falls off for large values of the scale around 1 GeV. We observe 
a cancellation between the scale dependence of the short- and long-distance contributions, 
i.e., the operator evolution in the quark picture is continued with the same pattern in the 
meson picture. The main uncertainty displayed in Fig. |3| originates from the dependence of 
the Wilson coefficients on A QCD . The uncertainty increases for very low values of the scale 
reflecting the poor perturbative behaviour expected at those scales especially for the large 
value of A QCD = 405 MeV. Within the (conservative) range of A QCD = 325 ± 80 MeV we 
considered, the value 405 MeV leads to the most distinct deviation from the experimental 
result which, however, does not exceed approximately 20 % of the observed value in the 
range 600 MeV < A c < 800 MeV where the minimum occurs and the dependence on the 
scale is weak. 

In Fig. |] we compare the results for Rea we obtain using the LO and NLO Wilson 
coefficients, respectively. In the HV scheme, for moderate values of A QCD introducing 
the NLO coefficients does not significantly affect the numerical values of the AI = 1/2 
amplitude which is found to be only slightly suppressed with respect to the LO result. 
The main effect of the NLO coefficients is that they further reduce the dependence on the 
matching scale. This statement does not hold within the NDR scheme. In this scheme, 
for A QCD = 245 MeV the effect of the NLO coefficients is also moderate but noticeably 
increases for large values of A QCD leading to a distinct suppression of the LO result. For 
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Figure 3: Rea (in units of MeV) with LO Zi for m s (l GeV) = 175 MeV and various values 
of A QCD as a function of the matching scale A c = ji. 

values of A QCD as large as 405 MeV both the HV and the NDR results rapidly diverge 
for low values of the matching scale (< 700 MeV) indicating the loss of perturbativity. 
Taking into account the fact that we do not incorporate the effects of higher resonances 
and cannot adopt too high values of the scale, a choice of A c around 700 - 800 MeV seems 
to be most appropriate. For A QCD = 325 MeV (245 MeV) the effect of the NLO coefficients 
is less pronounced, and scales as low as 600 -650 MeV (500 MeV), where the LO minimum 
occurs, appear to be acceptable. Above these scales the deviation of the NLO results 
from the experiment does not exceed 20 - 25 % of the experimental value. Moreover, the 
difference between LO and NLO (HV and NDR) values is moderate, of the order of at most 
20 - 25 % of the observed value.0 In all the cases the tendency for a large enhancement of 
10 The comparison of the LO and NLO coefficients should be used with caution as it partly originates 
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Figure 4: Rea (in units of MeV) with LO and NLO z; for m s (lGeV) = 175 MeV and 
various values of A QCD as a function of the matching scale A c = //. 

the required size remains present. 

In Fig. [5] we show the weak dependence of Reao (with LO Wilson coefficients) on the 
strange quark mass which arises from the matrix element of the gluon penguin operator 
[(Q6)o oc V m sl- We notice that the contribution from Q G to the AI = 1/2 amplitude 
for small values of the cutoff (~ 600 MeV) roughly varies between 10 - 20 % of the total 
value and significantly decreases for large values of A c . This behaviour is also found when 



from a change in the value of the QCD coupling for a chosen value of A^jg |li 
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Figure 5: Rea (in units of MeV) with LO Zi for A QCD = 325 MeV and various values of 
m s (l GeV) as a function of the matching scale A c = \i. 

the NLO coefficients are used. The effect of the remaining (penguin) operators is very 
small (below 1 % of the total result except for Q 4 which contributes at the level of —3 %). 
For comparison, in Fig. [^ we also show Reao calculated in the chiral limit. We observe 
that the result obtained in the chiral limit, for reasons explained above, is rather close 
to the numerically exact one, that is to say, the logarithmic and the finite terms in the 
non-factorizable corrections to the matrix elements are minor important provided that the 
matching scale is taken sufficiently large (A c > 500 -600 MeV). Finally, we note that the 
presence of the 770 does not affect the numerical values of the amplitudes (in the octet limit 
the numbers given in Table [8] change by less than 1 % ). 

In distinction to the A J = 1/2 amplitude, the AI = 3/2 amplitude depicted in Fig. || 
(with LO Wilson coefficients) is highly unstable. In addition, the numerical values lie well 
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below the measured value. The amplitude even changes sign [due to the large negative 
coefficient of the quadratic term in Eq. (ffB"|)]. The large uncertainty can be understood, 
as we already discussed above, from the fact that the two numerically leading terms, the 
tree level and the one-loop quadratically divergent term, have approximately the same size 
but opposite sign. On the one hand, this property is generally welcomed as it explains 
the origin of the suppression of the AI = 3/2 amplitude which turns out to be sufficiently 
suppressed whatever the particular chosen scale is between 600 MeV and 900 MeV. On the 
other hand, the large cancellation implies that the result will be significantly affected by 
higher order terms which are expected to be of the order of the one-loop logarithmic and 
finite terms. We note that the agreement with the experimental value is not improved 
in the chiral limit. We also notice that the numerical values depicted in Fig. |^ depend 
only weakly on the choice of A QCD . In Fig. |7| we compare the results for Rea2 we obtain 
using the LO and NLO Wilson coefficients, respectively. We observe that the effect of the 
NLO coefficients is negligible with respect to the large discrepancy between our results 
and the observed value. The small effect of the NLO coefficients indicates the validity of 
perturbation theory and further supports the supposition that the discrepancy is due the 
lack of accuracy in the low-energy part of the calculation. 

The typical size of higher order effects in the calculation of the hadronic matrix elements 
can be estimated in various ways. First, as we already mentioned above, one may replace 
in all NLO terms the coefficient l/F n by 1/Fk- The results obtained in this case [denoted 
by (b)] are shown in Figs. || and §. The AI = 1/2 amplitude is suppressed by approximately 
20 % with respect to the result we obtained using l/F n [denoted by (a)] and is even in 
better agreement with the observed value. The AI = 3/2 amplitude, on the other hand, 
is enhanced but still far too much suppressed. Another estimation of higher order effects 
can be done, as we explained above, by completely neglecting the imaginary part of the 
matrix elements (c). This suppresses Rea by a factor of cos<5q xp ~ 0.8 but does not 
affect Rea2- Similarly the absolute value of the amplitudes can be calculated by taking 
directly the imaginary part from Tables |3] and |] without using the experimental phases 
(d). This procedure suppresses Rea in the same way as in the previous case but largely 
re-stabilizes Rea2, indicating that the results obtained for the AI = 3/2 amplitude (unlike 
those obtained for Reao) indeed can be significantly affected by higher orders corrections. 
It is unlikely, however, that higher order terms alone can account for the large discrepancy 
between our result and experiment, and effects from higher resonances are also expected 
to be non- negligible for the small AI = 3/2 amplitude. Finally, the coefficient 1/F^ in the 
next-to-leading order terms can also be replaced by the bare coupling 1/ f as it was done in 
Ref . [pl| . Even though this would introduce an unphysical dependence on the factorizable 
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Figure 6: Rea 2 (in units of MeV) with LO Zi for various values of A QCD as a function of 
the matching scale A c = /x. 

scale, formally the difference also concerns higher order effects. Q We observe that this 
choice (e) leads to a result for Reao which is approximately scale independent. It also gives 
a more stable result for Rea 2 which, however, still is too much suppressed. 

In summary, in all cases we discussed above the AI = 1/2 amplitude is obtained 
around the measured value with an uncertainty of less than 25 % or in most cases even 
less than 15 %.Q The result for Reao is consequently solid and presumably could be 
significantly affected only by higher resonances. In view of the good agreement with the 

11 The relation between F^ and / is given in Eq. (62) of Ref. [EOf and we obtain / = 105, 112, 120, 
128, 136, 145 MeV for A c = 500, 600, 700, 800, 900, 1000 MeV, respectively. 

12 The only exception to this is the case where the large value of A qcd = 405 MeV is taken at LO or 
NLO (HV scheme) using a matching scale as high as ~ 1 GeV. In this (unfavourable) case the deviation 
from the observed value can be as large as 35-40%. 
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Figure 7: Rea2 (in units of MeV) with LO and NLO Zi for various values of A QCD as a 
function of the matching scale A c = //. 

experiment we obtained at the pseudoscalar level their effect a priori is expected to be 
small. The A/ = 3/2 amplitude, on the other hand, though showing the qualitatively 
correct behaviour of being suppressed with respect to the VSA result, emerges too much 
suppressed and is very unstable. However, higher order corrections to the matrix elements 
have been estimated large and could re-enhance it. In the same way higher resonances 
could easily enhance the result obtained at the pseudoscalar level. Vector mesons can be 
incorporated in a straightforward (however lenghty) way, and it would be very interesting 



27 



xlO 
0.8 

0.7 

0.6 

0.5 h 
0.4 



0.3 



0.2 



0.1 



- (a) 

- (b) 
• (c) 
-(d) 

■ (e) 

■ experiment 



_i i i i i i i i i_ 



_i i i i i i i_ 



600 650 700 750 800 850 900 950 1000 

A (MeV) 



Figure 8: Rea (in units of MeV) with LO z* for m,(lGeV) = 175 MeV and A QCD = 
325 MeV within different treatments of higher order corrections as explained in the text. 



to investigate their effect in the present calculation. This also would allow more safely to 
choose higher values for the matching scale for which the short-distance contributions are 
more reliable. 

We close this section by a brief review of several other attempts which have been 
made to explain the AI = 1/2 rule using different methods for the computation of the 
hadronic matrix elements. Interesting tendencies for an enhancement of the AI = 1/2 
channel were found in particular in Ref. JTTJ] by integrating out the quark fields in a 
gluonic background and in Ref. [|12 in the framework of QCD sum rules at the level 



of the inclusive two-point function. In Ref. |L3| quantitative results reproducing both the 
AI = 1/2 and AI = 3/2 channels were obtained adopting the point of view that in addition 
to 1/N C effects due to one-loop corrections (similar to those of Fig. ^) diquark states play 
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Figure 9: Rea 2 (in units of MeV) with LO Zi for A QCD 
ments of higher order corrections as explained in the text. 



325 MeV within different treat- 



an important role. The results for the AI = 1/2 amplitude obtained in the present 
approach suggest that there are no large diquark effects not already taken into account 
in the 1/N C corrections we calculated. The AI = 1/2 rule has also been investigated in 
the framework of chiral perturbation theory [I3| and the chiral quark model [|15| . At the 
present state of these methods the ratio 1/u = 22.2 cannot be predicted but is used to fit 
parameters of the models. Very recently the matrix elements relevant for the AI =1/2 
rule were studied in lattice QCD with improved statistics [ 16| . The authors used lowest- 
order chiral perturbation theory to relate the matrix elements (iT7i\Qi\K ) to (ix + \Qi\K + ) 
and (0\Qi\K°) calculated on the lattice. The ratio of the amplitudes computed in this 
way confirms the significant enhancement of the AI = 1/2 channel although systematic 
uncertainties preclude a definite answer. Whereas the AI = 1/2 amplitude is obtained 
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larger than the experimental value by approximately 40 % (quenched ensembleP], (3 = 6.0) 
the AI = 3/2 amplitude suffers from ambiguities in the choice of the meson mass due to 
the ignorance of higher order chiral corrections to the relation between Rea2 and the B^ 
parameter. Taking the meson mass M 2 = (m 2 K + m 2 )/2 and using the quenched value of 
Bk in the continuum limit the authors obtain a value for Rea2 which also over-estimates 
the data by approximately 40 %. The ratio of the amplitudes exhibits a strong dependence 
on the meson mass (see Fig. 11 of Ref. |ll|) due to the chiral behaviour of Rea2- In 



lattice perturbation theory unlike in analytical methods, the matching of the renormalized 
operators to the Wilson coefficients can be rigorously done, at least in principle (see e.g. 



Ref. [45| and references therein). On the other hand, analytical methods like the 1/N C 
approach followed in this paper allow for a direct evaluation of the K — > ttti amplitudes 
without the need of using reduction formulas to relate these amplitudes to the off-shell 
K — > re amplitudes (for this point see also Ref. |46| and references therein). 



While this paper was written an analysis of the AI = 1/2 rule was published |4"? 
which follows similar lines of thought as our work. In their analysis the authors used 
the 1/N C expansion in the chiral limit in the framework of chiral perturbation theory and 
the Extended Nambu-Jona-Lasinio model, respectively. We agree on the coefficients of 
the quadratically divergent terms in the 1/N C corrections to the matrix elements quoted 



therein. In the present analysis we did not investigate the method proposed in Ref. |47| to 
treat the scheme dependence appearing at the next-to-leading logarithmic order. 

5 K° -K° Mixing 

The contributions from short-distance physics to K° — K° mixing can be calculated from 
an effective AS = 2 hamiltonian, valid below the charm threshold, in which the heavy 



degrees of freedom are integrated out |32j , 

H?f f = 2 = F(m 2 t ,m 2 c ,M^V CKM )G F {a s (fi)]- 2/9 
where Oas=2 is the following four-quark operator: 



1 + ^T Js 



Oas=2 , (39) 



A S=2 = SL'fdL Sl^^l , (40) 

with a s (/x) being the QCD running coupling with three active flavors and J3 a renormal- 

ization scheme dependent coefficient appearing at the next-to-leading logarithmic order. 

13 Quantitative estimates of quenching effects on the coefficients of the chiral logarithms in the one- loop 
contributions to the K — > tt-k amplitudes were presented in Refs. J43], [44|. In Ref. J43[ finite volume effects 
on the lattice were also investigated. 
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.F(mf , m? c) M$y, V CKM ) is a known function of the heavy quark masses, the W boson mass, 
and CKM matrix elements. It incorporates the basic electroweak (box diagram) loop con- 
tributions [[S|, as well as, the perturbative QCD effects described through the correction 



factors T]i,r]2, 773 which have been calculated at the leading logarithmic ||, [49[ and the next 



to- leading logarithmic order [^, |33| . Terms depending on a s (/i) are factored out explicitly 
to exhibit the renormalization scale (and scheme) dependence of the coefficients which has 
to cancel the corresponding scale (and scheme) dependence of the hadronic matrix ele- 
ment of Oas=2 [0]- The short- distance hamiltonian for AS = 2 transitions in Eq. ( p9|) 



dominates the indirect CP violation in the neutral kaon system parameterized by e. Con- 
tributions to K° — K° mixing changing strangeness by two units through two AS = 1 
transitions at long distances which are relevant for the Kl — K$ mass difference |2S| are 
not considered in this article. 

The hadronic matrix element of 0&s=2 is usually parameterized in terms of the Bk 
parameter which quantifies the deviation from the value obtained in the vacuum saturation 
approximation: 

(K \O AS=2 (fi)\K°) = /^(/i)(/? o |0as=2|^°>vsa, (41) 

where 

(^ |Oas= 2 |^°)vsa = \F 2 K m\. (42) 



It is convenient to introduce the renormalization group invariant parameter [19], |50 
B K = B k (jjl) [a s (fi)]- 2/9 



1 + ^ J 3 

4tt 3 



1 (NDR) 



^3 = 7 , (43) 

1 (HV) 



in which the scale (and scheme) dependences of the long- and short- distance contribu- 
tions cancel within an exact realization of both perturbative and non-perturbative QCD. 
However, from the results for the AI = 3/2 K — ► nn amplitude discussed in the previous 
section we do not expect that the Bk we will obtain within the pseudoscalar approximation 
used in the low-energy calculation will exhibit a negligible dependence on the matching 
scale; the 27-plet operators which induce AS = 1 (AI = 3/2) and AS = 2 transitions 
are components of the same irreducible tensor under SU(3)l x SU(3)r, that is to say, 
to leading order in the chiral expansion the K° — K° amplitude can be related to the 
AI = 3/2 part of the K — > nn amplitude using 577(3) symmetry [51], [5^]. Consequently, 
we expect a similar pattern, i.e., a large negative quadratic term in the 1/N C corrections to 
the matrix element which partly cancels the tree level contribution and renders the result 
more sensitive to corrections from higher order terms and higher resonances. On the other 
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Figure 10: Factorizable contributions to the matrix element of the K° — K° mixing ampli- 
tude in the isospin limit. 



hand, we expect SU(3) breaking effects in AS* = 2 transitions to be more pronounced than 
in AS = 1 transitions |k|. In the following we will see that the 1/N C expansion restricted 
to the pseudoscalar mesons indeed leads to a significantly scale dependent result for B K . 
However, the scale dependence is less pronounced than the one of the AI = 3/2 amplitude 
due to corrections beyond the chiral limit. Finally, as we already discussed above, the 
low-energy calculation does not allow any reference to the renormalization scheme depen- 
dence. Nevertheless, a comparison of the Bk parameter obtained from the LO and NLO 
coefficient function of 0as=2 can be used to to test the validity of perturbation theory and 
to estimate the uncertainties arising from the short-distance part. 



5.1 Factorizable Loop Corrections 

To obtain the factorizable non-perturbative corrections to the AS = 2 transition we have 
to calculate the diagrams in Fig. [Hi Using the chiral representation of the quark current 



in Eq. ([Tq) and reducing the result to the basic integrals listed in Appendix B of Ref. |20 
we obtain the unrenormalized (bare) matrix element: 



(K \O, 



A '7=2 1 -f^ )( ) 



2 f2 



m K f 



t . 1QL 5 2 

1 H TT^rn 



P 



K 



9/ 2 



(a + 2 bf h [mj +2 (a-bf h [m n >] + 18 h [m K ] + 9 h [m n 



(44) 



with a and b defined in Eq. (|i~5l) . Multiplying Eq. (f4"4"D with Z K , i.e., including a factor 



7-1/2 
J K 



for each external kaon field (compare Eqs. (16) and (59) of Ref. 



we arrive at 



(K \O AS=2 \KY = m\f 



1 H — —m 



P 



K 



up 



9 h [m n ] + 18 h [m K ] + (a + 2 bf h [m v ] + 2 (a-bf h [m rf 



(45) 



Comparing Eq. ( |4"5D with Eqs. (26) and (63) of Ref. [^0| we observe that the correction 
factor in the brackets which is due to the higher order (factorizable) contributions to the 
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Figure 11: Non-factorizable contributions to the matrix element of the K° — K° mixing 
amplitude in the isospin limit. 

matrix element is completely absorbed (including the finite terms) in the renormalization 
of the kaon decay constant, as it is required by current conservation, leading to the final 
result for the (renormalized) factorizable matrix element 

{K°\0^ S =2\K%. ) = m\F 2 K . (46) 

Eq. ( fl6|) represents the large- N c limit for the K° — K° matrix element, i.e., B K C ~*°° = 3/4, 
to be compared with the VSA value one. 

5.2 Non-factorizable Loop Corrections 

The 1/N C corrections to Eq. (pj|) can be calculated from the non-factorizable loop diagrams 
depicted in Fig. [11]. We determine the loop momenta along the lines developed in Section 2, 
that is to say, by associating the cutoff to the effective color singlet boson connecting the 
two currents. The simple structure of the non-factorizable diagrams makes it possible to 
specify the complete analytic result for the matrix element in terms of loop integrals. In 
the SU(2) limit the expression in which the integrals are reduced to the basic ones reads 

(K \O AS=2 \K°) NF = ^ + l^m 2 K -2p 2 K -(x2 + X3))h[m K } 

- t:(X2 + Xs + 2rn% + 2p 2 K )m 2 K I 3 [m K , m K , 0] - -{p 2 K + m^)J 2 [m 7r ,p ft :] 

13 3 

+ -h[m n ,p K ] + -cos 2 6 I 4 [m v ,p K ] + - sin 2 6 1 4 [my , p K \ . (47) 

Here we replaced a and b by the 77 — 77' mixing angle 9 and explicitly distinguished be- 
tween the masses coming from the external kaon momentum, the explicit mass term in 
the lagrangian, and the propagators in the loops. In addition to the logarithmically and 



33 



quadratically divergent integrals (/]_, ^h) listed in Appendix B of Ref. |20| Eq. fl47|) con- 
tains the integral I4 which exhibits a quartic dependence on the cutoff. Following the steps 



discussed in Ref. |K| we can give the analytic expression for J 4 in terms of a Taylor-series: 

,2 



h[m,p] 



(2vr) 4 

1 

16^2 



d 4 q 



1 



Q 



(q — p) 2 — m 2 



' 2 A ' + m * 



A„ — m lo£ 



m 2 



+ 



2 2 
p m 



(A c + m 2 ) 2 
p 4 A° 



^ + A 2 m 2 - (A 2 + m 2 ) 2 log (l + ^f 



6(A 2 + m 2 ) 4 



(A 2 - 2m 2 ) + 



p 6 A^m 



2(A 2 + m 2 ) 



A: 



—m 
3 



+ 0(p s ) . (48) 



We note that the logarithmically divergent integral I3 in Eq. ( [4"TD only appears with van- 
ishing external momentum and therefore can be largely simplified compared to the general 
expression in Eq. (75) of Ref. [20 1 . From Eq. (47) one can easily calculate the divergent 
terms. Taking the external momentum on-shell we obtain 



(K \O AS=2 \K ) 



0\NF 



m\Fl 



3A 2 



(4tt) 2 F1 



(4m t- — 2m 2 Y m 2 + mt) , . 9 

+ - K lA * \ n> log A 2 + ■ 

{^) 2 F 2 K m 2 K 



(49) 

where the tree level result is factored out and the ellipses denote the finite terms we do not 
specify analytically. We observe that the quartic dependence on the cutoff is cancelled as 
required by chiral symmetry. 

To illustrate the effect of the modified momentum routing we also recalculate the non- 
factorizable loop contributions in the approach used by Bardeen et al. [EB| who associated 
the cutoff to the momentum of the virtual meson in the loop diagrams (see also the dis- 
cussion in Ref. [p0f ): 

1 
12 



{K \U AS=2 \K ) BBC 



2 (x2 + X3 - 2rn 2 K ) h[m K ] 
+3 (m 2 K + m 2 ) I x \rn v ] + 9 cos 2 6 {m\ + m 2 ) I x [mj + 9 sin 2 
x {m\ + m 2 ,) I^] + 2m 2 , ( %2 + Xs + 4m 2 ,) I 3 [m K , m K , 0] 



where the external momentum is already taken on-shell. For comparison with Eq. ([47|) in 
Eq. flBTf ) we included the small effect of the singlet rj . Solving the integrals we obtain the 



(50) 
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A c 


0.5 GeV 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


(@AS=2) tree 
(Oas=2) a 2 
(C?AS=2)log + fin 


3.14 

-1.17 

0.57 


3.14 

-1.68 

0.76 


3.14 

-2.29 

0.96 


3.14 

-2.99 

1.15 


3.14 

-3.78 

1.32 


3.14 

-4.67 

1.49 


(0 AS =2) 


2.54 


2.22 


1.81 


1.30 


0.68 


-0.04 


Bk(Ac) 


0.61 


0.53 


0.43 


0.31 


0.16 


-0.01 



Table 9: Different contributions to the hadronic matrix element of Oas=2 (in units of 
10 9 • MeV 4 ) and Bk, shown for various values of the cutoff A c . 



divergent part of the non-factorizable loop corrections 

(K°\Oas=2\K°) 



-0\NF 

'BBG 



m K F K 



2kl Um% - 2mtmi + mi) , . , 

+ K ,. *o o " ; log A* + 



(4tt) 2 F1 



K 



w 2 n< 



K 



(51) 

to be compared with Eq. (|4*9|). We note that the results obtained in both calculations differ 
with respect to the quadratic cutoff dependence, as well as, with respect to the finite terms 
we do not give explicitly here for brevity. 



5.3 Numerical Results 



As a numerical input we use the values listed in Section P . In Table [] we show our results 
for the K° — K° matrix element and Bk(A. c ) obtained in the full calculation, i.e., including 
the effect of the r]o in Eq. fl4~7|) . In Fig. [12| we depict the renormalization group invariant 
parameter Bk calculated with the leading order Wilson coefficient. 

The decrease of Bk(A c ) with A c = /i is qualitatively consistent with the /i dependence 
of the coefficient function in Eq. (|43|) , that is to say, the long-distance evolution counteracts 
the evolution in the short- distance domain. This property is due to the presence of the 
quadratic terms in the 1/N C corrections which compensate for the (weaker) increase of 
the logarithmic terms. However, the decrease is found to be significant, and the scale 
dependence largely exceeds what is required to have an exact cancellation of both evolutions 
over a large range of the scale. As a result an acceptable stability of B^ is obtained only 
for low values of A c ~ 500 - 600 MeV. The small values of Bk depicted in Fig. ^ (lower 
set of curves) come from the negative coefficient of the quadratic term in Eq. (|4"9| ) which is 
found to be enhanced by a factor of 3/2 compared to the result of Ref. p6j. This coefficient 
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Figure 12: Br with LO Wilson coefficient for various values of A QCD as a function of the 
matching scale A c = fi. The lower set of curves shows the results of the present analysis, 
the upper set allows a comparison with Ref. [EB . 



is the same as the one of the AI = 3/2 K —* titt amplitude except for ST/ (3) breaking 
effects (responsible for Fk ^ F„) which reduce the negative slope of Bk- As can be seen 
from Table |], the difference between the exact result and the one obtained in the chiral 
limit (i.e., in the absence of chiral logarithms and finite terms) is more pronounced than 
in the case of the K — > tttt amplitudes. This is due mainly to the numerical coefficient of 
the leading term (~ frig) in front of the logarithm in Eq. (|49|) which as expected is found 
larger in AS" = 2 transitions than in AS = 1 transitions. Because of the large positive 
coefficient the logarithmic term re-stabilizes Bk sizably with respect to the result obtained 
in the chiral limit. This also explains why the Bk parameter even if significantly scale 
dependent is much more stable than the AI = 3/2 amplitude. The finite terms beyond the 



3(3 



logarithms in Eq. (47) [i.e., beyond the log(l + A;?/m 2 ) terms] give a negative contribution 



to -Brt(A c ) roughly between —0.05 and —0.08 for A c around 600 -900 MeV. Consequently, 
they are non-negligible in particular for large values of the scale where the cancellation 
between the tree level and the quadratic terms is large. Finally, we note that the presence 
of the f]o does not significantly affect the numerical values of the K° — K° matrix element 
(in the octet limit the numbers given in Table || change by less than 3 % ). 

To illustrate the effect of the momentum routing, in Fig. we also show Bk obtained 
from Eq. ( pip (upper set of curves). We use the same set of parameters as in Table |9| and 
also include the rjo. Comparing the two results we notice that 5^(A C ) calculated within the 
modified momentum routing lies below the values found in the previous approach. Match- 
ing the long-distance results with the short- distance contribution we observe that the Bk 
parameter obtained in the present analysis exhibits a significantly stronger dependence on 
the matching scale. However, as we already discussed above, the quadratically divergent 
terms (and the finite terms) depend on the way we define the integration variable inside the 
loop. This can be seen from the different numerical factors in front of the quadratic terms 



in Eqs. (fE|) and ([51]). Therefore we are forced to find a direct link between the short- and 
long-distance part of the calculation, as it is done by keeping track of the effective color 
singlet boson in both parts of the calculation. A consistent matching is then obtained by 
assigning the same momentum to the color singlet boson at long and short distances and by 
identifying this momentum with the loop integration variable (see Section |2|). This prop- 
erty is absent in the previous approach. The modification unambiguously determines the 
coefficient in front of the (quadratically and logarithmically) divergent terms and allows us 
to identify the ultraviolet cutoff of the long-distance terms with the short-distance renor- 
malization scale fi. Therefore we advocate the use of the modified matching prescription, 
even though the stability of our result is rather poor. The satisfactory stability obtained in 
Ref. [Q on the other hand is somehow inconclusive, as there is no underlying argumenta- 
tion determining the quadratic terms. Our result also implies that the uncertainties due to 
the idealized identification of the cutoff A c with the upper limit of the meson momentum 
in the loop in Ref. [T2^] might have been underestimated. In a complete meson theory 
the dependence on the momentum routing should be absent. However, as long as we are 
working in an effective low-energy approach as chiral perturbation theory we have to pay 
attention to this point. 

Numerically, we find a range of acceptable stability in the energy regime from 500 MeV 
to 700 MeV (see Fig. [E|) leading to values for Bk in the range of 0.4 < Bk < 0.6. The lower 
bound corresponds to a value of A QCD = 405 MeV, whereas the upper bound corresponds 
to A QCD = 245 MeV. Comparing our result with the one of Ref. p6| we observe a tendency 
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Figure 13: B K with LO and NLO Wilson coefficient for various values of A QCD = A^ as 
a function of the matching scale A c = /i. For each value of A QCD the lower (intermediate, 
upper) curve shows the LO (HV, NDR) result. 

for Bk to be decreased to values below 0.6. This behaviour is due to the enhancement 
of the negative coefficient in front of the quadratic term in the 1/N C corrections to the 
K° — K° matrix element and, to a smaller extend, also due to the finite terms omitted 
in Ref. [J26| . However, our result suffers from a sizable dependence on the matching scale 
which precludes a precise answer. 

In Fig. [H| we compare the results for Bk we obtain with the LO and NLO coefficient 
function. For A QCD = 325 MeV in the HV scheme, introducing the NLO coefficient does 
not significantly affect the numerical values of the Bk parameter which is found to be only 
slightly enhanced with respect to the LO result. In the NDR scheme, the effect of the NLO 
coefficient is also moderate for large values of the scale but noticeably increases for low 
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values. For very low values of A c ~ 500 MeV the NLO result can differ from the LO one 
by as much as 0.2. However, for these scales the scheme dependence increases rapidly and 
it is desirable to take (at least) a matching scale around 600 -650 MeV where Bk is still 
relatively smooth and roughly varies between 0.45 and 0.6. For A QCD = 245 MeV in both 
the HV and NDR schemes a matching scale as low as 500 MeV appears to be acceptable, 
and within the range A c ~ 500 -650 MeV Bk is obtained between 0.5 and 0.7. On the 
other hand we observe that the pseudoscalar approximation would simply fail if A QCD was 
found as large as 405 MeV, as a satisfactory perturbative behaviour is obtained only for 
A c > 700 MeV, that is to say, for values of the scale where the stability of Bk is found 
to be poor. 

In summary, for values of A QCD > 350 MeV an estimate of Bk is hindered by the loss of 
perturbativity in the range where the pseudoscalar approximation is expected to be valid, 
and for lower values of A QCD (taking into account the scheme dependence) our calculation 
favours low values of Bk in the range 

0.4 < B K < 0.7. (52) 

However, a satisfactory smooth behaviour is obtained only in a narrow range of the cutoff 
and, in addition, for values of the cutoff as low as the kaon mass or just above. Therefore 
the incorporation of higher resonances is clearly required as for the AI = 3/2 K — ► tttt 
amplitude discussed above. On this issue, the analysis of the Bk parameter is similar to 
the one of the AI = 3/2 amplitude, even if numerically the matching obtained for Bk is 
better than the one obtained for the AI = 3/2 amplitude. 

The K° — K° system has been studied in the past with various methods leading to 
different results for B K - The present status of quenched lattice calculations [54-57] has 
been reviewed in Ref. |58|. The value reported by the author is B K = 0.86 ± 0.06 ± 0.06. 



Very recently the JLQCD Collaboration has presented a new analysis based on chiral 
Ward identities to non-perturbatively determine the mixing coefficients of the AS = 2 



operator ||59| . The numerical results given in Ref. [[59] are in agreement with the lattice 
calculations quoted above. In the chiral quark model a value as high as Bk = 1.1 ± 0.2 
has been obtained ||15|| . Lower values for Bk have been found in the QCD hadronic duality 



approach @ (B K = 0.39 ± 0.10), by using ST/ (3) symmetry and PCAC [H (~ 1/3), 



or using chiral perturbation theory at next-to-leading order |H[] (0.42 ± 0.06). QCD sum 



rules give results around B K = 0.5-0.6 with errors in the range of 0.2-0.3 |62|, |63|. One 
might note that a value for Bk significantly below 0.7 requires simultaneously high values 
of |Kift/Vd,| and \Vd>\ to be able to fit the experimental value of e |I9| . Finally, we note 



that the Bk parameter was also investigated in the framework of the 1/N C expansion in 
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Ref. 1 50fl . In this work the matching was not performed at the level of the K° — K° matrix 



element but at the level of a related 2-point Green function. Numerically, the matching was 
found unsatisfactory good. We agree with this conclusion, as we discussed above, although 
in Ref. ]50| the quadratic dependence on the UV cutoff was obtained in disagreement 



with the present analysis due to the use of a different momentum routing. This has been 



corrected very recently in Ref. [f|7]], and we agree with the results for the 1/N C corrections 
to the K° — K° matrix element obtained there in the chiral limit. In the present paper we 
investigated the corrections beyond the chiral limit and found that they are sizable. On the 
other hand, the authors of Ref. ||7|] investigated higher order corrections calculated in the 
framework of the Extended Nambu- Jona-Lasinio model. As a result they obtained a better 
stability of the Bx parameter. This shows that corrections from higher order terms and 
higher resonances are expected to be large. Nevertheless the values of Bk we obtained in 
this analysis by performing a full calculation at the pseudoscalar level are meaningful and 
can be considered as reference values for further investigations incorporating the effects of 
higher resonances. 

6 Conclusions 

The 1/N C approach developed in Refs. fl8 , p6| when modified along the lines of Ref. p0 | 



leads to interesting results in the current- current sector of the AS = 1 and in the AS = 2 
transitions. The main result of the present analysis is an additional enhancement of the 
AI = 1/2 channel in the K — ► tttt amplitudes. This channel has been found sufficiently 
enhanced, in good agreement (with an accuracy of 80 to approximately 100 %) with the 
experiment, and widely stable over a large range of values of the matching scale roughly 
between 600 MeV and 900 MeV. It is certainly premature to say that the dynamical mech- 
anism behind the AI = 1/2 enhancement is completely understood. An agreement at the 
level obtained in the present analysis a priori is not expected in an effective theory with 
only pseudoscalar mesons taken into account. Nevertheless we believe that the additional 
enhancement reported here is a further important indication that the 1/N C approach can 
account for the bulk of the AI = 1/2 amplitude. This statement is also supported by 
the fact that higher order corrections both of short- distance origin and of long-distance 
origin at the pseudoscalar level, as we discussed above, are not expected to largely affect 
the size of the AI = 1/2 enhancement. The agreement with the experiment also tends to 
show that the origin of the long-distance enhancement has to be found at the level of the 
pseudoscalar mesons and at energies below the rho mass or even below the kaon mass. Cer- 
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tainly this has to be checked explicitly incorporating at least the effects of vector mesons. 
We also believe that the 1/N C approach can account for the bulk of the suppression of the 
A/ = 3/2 channel. For this channel, however, the approximations made in the present 
analysis fell short of the desired accuracy. In particular, a large scale dependence has been 
found clearly requiring the incorporation of higher order terms and/or higher resonances. 
We note that the scale behaviour of the ratio of the two isospin amplitudes is dominated 
by the one of the A/ = 3/2 channel, and therefore it leads to a comparable uncertainty. 
Similarly, the Bk parameter suffers from a sizable dependence on the matching scale. Our 
calculation favours very low values of the scale (< 700 MeV) leading to values for Bk in 
the range of 0.4 < Bk < 0.7. However, the large uncertainties associated with this result 
preclude a definite answer, and also make the incorporation of higher order terms and 
higher resonances very desirable. 



Acknowledgments 

We wish to thank Emmanuel Paschos for many valuable discussions. We also thank William 
Bardeen, Johan Bijnens, Jorge Fatelo, Jean-Marc Gerard, and Gino Isidori for helpful 
comments. We are very thankful to Matthias Jamin for providing us with the numerical 
values of the Wilson coefficients used in Section 4 of this article. This work was supported 
in part by the Bundesministerium fur Bildung, Wissenschaft, Forschung und Technologie 
(BMBF), 057D093P(7), Bonn, FRG, and DFG Antrag PA-10-1. 



41 



A Numerical Values of the Wilson Coefficients 



In this appendix we list the numerical values of the LO and NLO (HV and NDR) Wilson 
coefficients for AS = 1 transitions used in Section [L2"[ These values were supplied to us by 
M. Jamin. Following the lines of Ref. || the coefficients Zi are given for a 10-dimensional 
operator basis {Qi, . . . , Qw}. Below the charm threshold the set of operators reduces to 
seven linearly independent operators [see Eqs. (§)-(0)] with 



Qi 



-Qi + Q2 + Qi 



Q9 = ~Qi-~Q 3 , Qio = \Qi + Q2-\Qz- (53) 



At next-to-leading logarithmic order in (renormalization group improved) perturbation 
theory the relations in Eq. (|53|) receive 0(a s ) and 0(a) corrections 0, |i~9j . In the present 
analysis we use the linear dependence at the level of the matrix elements (Qi)i, i.e., at the 
level of the pseudoscalar representation where modifications to the relations in Eq. ( |53"D 
are absent. We note that the effect of the different treatment of the operator relations at 
next-to-leading logarithmic order which is due to the fact that in the long-distance part 
there is no (perturbative) counting in a s is numerically negligible. 

The following parameters are used for the calculation of the Wilson coefficients: 



M w = 80.2 GeV, 



sin 2 6 



w 



0.23, 



a 



1/129, 



m t 



170 GeV, 



m b [m b ) 



4.4 GeV, 



m c {m c 



1.3 GeV. 



M 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-0.937 


-0.826 


-0.748 


-0.690 


-0.645 


Z2 


1.576 


1.491 


1.433 


1.391 


1.359 


Z3 


0.016 


0.011 


0.007 


0.005 


0.003 


z 4 


-0.037 


-0.027 


-0.019 


-0.014 


-0.009 


Z5 


0.011 


0.008 


0.006 


0.004 


0.003 


Z6 


-0.045 


-0.031 


-0.021 


-0.015 


-0.010 


z 7 /a 


0.023 


0.017 


0.012 


0.008 


0.005 


z$/a 


0.007 


0.004 


0.002 


0.001 


0.0004 


z 9 /a 


0.027 


0.019 


0.013 


0.009 


0.006 


z w /a 


-0.006 


-0.003 


-0.002 


-0.001 


-0.0004 



Table 10: AS = 1 LO Wilson coefficients for A 



QCD 



245 MeV. 



42 



IX 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-1.192 


-1.010 


-0.893 


-0.811 


-0.748 


Z2 


1.779 


1.632 


1.541 


1.479 


1.433 


Z 3 


0.025 


0.016 


0.010 


0.007 


0.004 


Za 


-0.054 


-0.036 


-0.026 


-0.018 


-0.012 


z 5 


0.015 


0.011 


0.008 


0.006 


0.004 


z 6 


-0.070 


-0.044 


-0.029 


-0.019 


-0.013 


z 7 /a 


0.033 


0.023 


0.017 


0.012 


0.008 


z$/a 


0.012 


0.006 


0.003 


0.001 


0.001 


z 9 /a 


0.040 


0.027 


0.019 


0.013 


0.008 


zio/a 


-0.010 


-0.005 


-0.003 


-0.001 


-0.001 



Table 11: AS = 1 LO Wilson coefficients for A 



QC'D 



325 MeV. 



fi 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-1.576 


-1.246 


-1.065 


-0.947 


-0.861 


Zl 


2.104 


1.824 


1.676 


1.582 


1.517 


z 3 


0.041 


0.023 


0.014 


0.009 


0.006 


za 


-0.082 


-0.051 


-0.034 


-0.023 


-0.015 


z 5 


0.022 


0.015 


0.010 


0.007 


0.005 


Z& 


-0.119 


-0.066 


-0.041 


-0.026 


-0.016 


z-j/a 


0.044 


0.031 


0.022 


0.015 


0.010 


z 8 /a 


0.024 


0.010 


0.005 


0.002 


0.001 


z 9 /a 


0.056 


0.037 


0.025 


0.017 


0.011 


zio/a 


-0.017 


-0.008 


-0.004 


-0.002 


-0.001 



Table 12: AS = 1 LO Wilson coefficients for A 



QC'D 



405 MeV. 
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IX 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-0.668 


-0.578 


-0.516 


-0.470 


-0.435 


Z2 


1.391 


1.326 


1.282 


1.252 


1.229 


Z 3 


0.038 


0.023 


0.016 


0.012 


0.009 


Za 


-0.088 


-0.059 


-0.043 


-0.032 


-0.025 


z 5 


0.007 


0.009 


0.008 


0.007 


0.006 


z 6 


-0.102 


-0.064 


-0.044 


-0.032 


-0.025 


z 7 /a 


0.018 


0.012 


0.008 


0.006 


0.005 


z$/a 


0.069 


0.039 


0.024 


0.015 


0.009 


z 9 /a 


0.045 


0.029 


0.020 


0.014 


0.010 


zio/a 


-0.032 


-0.021 


-0.014 


-0.009 


-0.006 



Table 13: AS = 1 NLO Wilson coefficients (NDR) for A 



QC'D 



A. 



(-1) 



245 MeV. 



n 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-0.898 


-0.739 


-0.644 


-0.579 


-0.531 


Zl 


1.569 


1.444 


1.373 


1.326 


1.292 


z 3 


0.033 


0.019 


0.012 


0.007 


0.005 


za 


-0.060 


-0.038 


-0.025 


-0.017 


-0.011 


z 5 


0.012 


0.008 


0.006 


0.004 


0.003 


Z& 


-0.060 


-0.036 


-0.024 


-0.016 


-0.010 


z-j/a 


-0.005 


-0.005 


-0.004 


-0.004 


-0.003 


z 8 /a 


0.046 


0.027 


0.017 


0.011 


0.007 


z 9 /a 


0.023 


0.012 


0.006 


0.003 


0.001 


zio/a 


-0.038 


-0.024 


-0.016 


-0.010 


-0.007 



Table 14: AS = 1 NLO Wilson coefficients (HV) for A 



QCD 



A^ 

MS 



245 MeV. 
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IX 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-0.805 


-0.712 


-0.623 


-0.558 


-0.509 


Z2 


1.495 


1.424 


1.359 


1.312 


1.278 


Z 3 


0.095 


0.046 


0.027 


0.018 


0.013 


Za 


-0.193 


-0.104 


-0.068 


-0.048 


-0.035 


z 5 


-0.019 


0.005 


0.009 


0.009 


0.008 


z 6 


-0.261 


-0.121 


-0.072 


-0.049 


-0.035 


z 7 /a 


0.039 


0.025 


0.018 


0.014 


0.011 


z$/a 


0.181 


0.079 


0.042 


0.024 


0.014 


z 9 /a 


0.086 


0.054 


0.036 


0.025 


0.018 


zio/a 


-0.056 


-0.034 


-0.021 


-0.013 


-0.008 



Table 15: AS = 1 NLO Wilson coefficients (NDR) for A 



QC'D 



A. 



(-1) 



325 MeV. 



n 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-1.381 


-1.011 


-0.827 


-0.716 


-0.640 


Zl 


1.982 


1.662 


1.513 


1.427 


1.370 


z 3 


0.090 


0.040 


0.022 


0.013 


0.007 


za 


-0.129 


-0.068 


-0.041 


-0.026 


-0.016 


z 5 


0.016 


0.011 


0.008 


0.006 


0.004 


Z& 


-0.137 


-0.067 


-0.038 


-0.024 


-0.014 


z-j/a 


-0.008 


-0.003 


-0.002 


-0.002 


-0.002 


z 8 /a 


0.107 


0.050 


0.027 


0.016 


0.010 


z 9 /a 


0.052 


0.027 


0.015 


0.009 


0.005 


zio/a 


-0.077 


-0.042 


-0.025 


-0.016 


-0.010 



Table 16: AS = 1 NLO Wilson coefficients (HV) for A 



QCD 



A^ 

MS 



325 MeV. 
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/! 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-0.176 


-0.795 


-0.738 


-0.657 


-0.592 


Z2 


0.911 


1.485 


1.444 


1.384 


1.336 


Z 3 


0.350 


0.108 


0.052 


0.030 


0.019 


Za 


-0.637 


-0.218 


-0.117 


-0.074 


-0.050 


z 5 


-0.318 


-0.027 


0.004 


0.009 


0.009 


z 6 


-1.172 


-0.288 


-0.132 


-0.077 


-0.050 


z 7 /a 


0.119 


0.042 


0.029 


0.023 


0.018 


z 8 /a 


0.699 


0.185 


0.081 


0.042 


0.023 


z 9 /a 


0.132 


0.089 


0.059 


0.040 


0.029 


zio/a 


-0.077 


-0.054 


-0.033 


-0.020 


-0.012 



Table 17: AS = 1 NLO Wilson coefficients (NDR) for A 



QC'D 



A. 



(•i) 



405 MeV. 



n 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Z\ 


-2.603 


-1.494 


-1.102 


-0.901 


-0.778 


Zl 


3.138 


2.084 


1.739 


1.573 


1.475 


z 3 


0.370 


0.102 


0.044 


0.023 


0.012 


za 


-0.403 


-0.140 


-0.072 


-0.042 


-0.025 


z 5 


0.035 


0.014 


0.010 


0.007 


0.005 


Z& 


-0.463 


-0.141 


-0.067 


-0.037 


-0.021 


z-j/a 


-0.063 


-0.009 


-0.002 


-0.001 


-0.001 


z 8 /a 


0.342 


0.105 


0.048 


0.026 


0.014 


z 9 /a 


0.111 


0.051 


0.028 


0.016 


0.009 


zio/a 


-0.179 


-0.078 


-0.042 


-0.024 


-0.014 



Table 18: AS = 1 NLO Wilson coefficients (HV) for A 



QCD 



A^ 

MS 



405 MeV. 
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